Noname manuscript No. 

(will be inserted by the editor) 



Molecular data and radiative transfer tools for ALMA 



Floris van der Tak • Michiel Hogerheijde 



o 
o 



Received: 3 1 January 2007 / Accepted: 



> 

in 
oo 
m 
(N 
O , 
l> . 

o 

Oh 
I 

o 



Abstract This paper presents an overview of several mod- 
eling tools for analyzing molecular line observations at sub- 
millimeter wavelengths. These tools are already proving to 
be very valuable for the interpretation of data from current 
telescopes, and will be indispensable for data obtained with 
ALMA. The tools are: (1) the Leiden Atomic and Molecular 
DAtabase (LAMDA), a collection of spectroscopic data and 
collisional excitation rates; (2) RADEX, an on-line and off- 
line program to calculate non-LTE excitation and emission 
from a homogeneous medium, based on the escape proba- 
bility approximation; (3) RATRAN, an accelerated Monte 
Carlo program to solve molecular excitation and radiative 
transfer in spherical and cylindrical symmetry. The paper 
presents examples of how to use these tools in conjunction 
with existing data reduction packages to quantitatively inter- 
pret submillimeter single-dish and interferometric observa- 
tions. The described tools are publically available at 

http : //www . strw . leidenuniv.nl/ ^moldata. 

The paper concludes with a discussion of future needs in 
the fields of molecular data and radiative transfer 



1 Introduction 

Observations of spectral lines at radio, (sub)millimeter and 
infrared wavelengths are powerful tools to investigate the 
physical and chemical conditions in the dilute gas of astro- 
nomical source s where t hermodynamic equilibrium is a poor 
approximation ( Genze]|fl99 1 ; Black 2000). To extract astro- 
physical parameters from the data, the excitation and opti- 
cal depth of the lines need to be estimated, for which var- 
ious metho ds may be used, depending on the available ob- 
serva tions ( Van Dishoeck & Hogerheiidell999l ; IVan der Tak 
2005h. 



With the advent of ALMA, the need for efficient molecu- 
lar line modeling tools will become more pressing than ever. 
This paper presents several such tools which we have devel- 
oped and which we expect will become quite important for 
the analysis of ALMA data. 
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2 LAMDA 



The L eiden Atomic and Molecular Database (LAMDA) dSchoier et al, 
2005h contains spectroscopic and collisional data for 24 species 
of astrophysical interest. Most of these species are molecules, 
but data files for the fine structure lines of C, O and C+ are 
also provided. The data files contain energy levels, statisti- 
cal weights, Einstein A— coefficients and collisional rate co- 
efficients. Available collisional data from quantum chemi- 
cal calculations and experiments have been extrapolated to 
higher energies (up to E/k ^ lOOOK). The format of the 
data files is such that they can be read directly into RADEX 
and RATRAN. 
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3 RADEX 



Description RADEX (IVan der Tak et al .I2OO7I) IS a computer 
program to calculate the intensities of atomic and molecu- 
lar lines produced in a uniform medium, based on statisti- 
cal equilibrium calculations involving collisional and radia- 
tive processes. The treatment includes radiation from back- 
ground sources; optical depth effects are calculated using 
the escape probability approximation. The user has a choice 
of geometries: uniform sphere, expanding sphere, or plane- 
parallel slab. The expanding sphere solution is analogous to 
the well-known LVG method. 

Availability The program is available both as an on-line cal- 
culator for quick estimates and as an off-line package for 
detailed analysis of multi-line datasets. 

Typical Uses ( 1 ) Multi-line observations of a region are avail- 
able, and an estimate is needed of the average temperature, 
density, and column density. (2) Based on estimates of the 
typical conditions in a source, estimates are made of the ex- 
pected line intensities, for example to prepare an observing 
proposal. 

Example Ratios of submillimeter lines of HCN and HCO+ 
trace densities between 10^ and 10** cm^^^ and temperatures 
between 20 and 200 K. By observing 3^ different transi- 
tions, density and temperature can be found (Fig.[ril. 



4 RATRAN 



Description RATRAN ("Ho gerheiide & van der Takl2000l) is 
a computer program to simulate non-LTE line radiative trans- 
fer in spherical or cylindrical symmetry. The program uses 
the accelerated Monte Carlo method to represent the radia- 
tion field, which offers both speed and reliability. 

Within these geometries, the source model can have any 
density, temperature, abundance, and velocity distribution. 
The program produces synthetic spectral line cubes for a 
given source distance and orientation. Third-party software 
such as Miriad or IDL readily produce higher-level output 
such as beam-convolved line fluxes or even simulated inter- 
ferometric visibilities. 

Availability The spherically symmetric version of the pro- 
gram is available as an off-line package. The cylindrically 
symmetric version is available on a collaborative basis. 

Typical Use The user has obtained molecular line observa- 
tions of an object, and wishes to quantitatively test a par- 
ticular model. With RATRAN, he produces synthetic data 
and compares these to the observations by eye, or uses a 
minimization to arrive at an optimum fit. 



Example iBrinch et alj (l2007b compare observations of the 
Young Stellar Object L 1489 IRS to a model with infalling 
and rotational motions. The data set consists of millime- 
ter interferometric observations, near-infrared scattered light 
imaging, and CO ro-vibrational lines. Comparison of data 
and model uses a Voronoi optimisation algorithm. 



5 Future molecular data needs 

Of the 1 34 molecules which are known in space, collisional 
rate coefficients (CRCs) exist only for twenty-four (§|2]i. Ac- 
curate CRCs are important for reliable estimates of molecu- 
lar column densities, and essential to estimate other param- 
eters such as kinetic temperature and volume density. Re- 
cently several important molecules for which only inaccu- 
rate CRCs were known, based on outdated, low-quality po- 
tential energy surfac es, have been re- calculated, such as CS, 
HC3N, and CH3OH (lDubernetl2005l) . 

Calculation of rate coefficients is a highly specialized 
and time consuming effort. New results appear in the liter- 
ature at a rate of a few molecules per year, which is at or 
below the detection rate of new interstellar molecules. The 
point where CRCs are known for all interstellar molecules 
may thus never be reached. In order to make the best pos- 
sible use of scarce human and computational resources, we 
have selected those molecules for which the calculation of 
new CRCs would have the most astrophysical impact. 

The main criterion to assign high priority to a molecule 
is its importance in constraining astrophysical or astrochem- 
ical parameters beyond its own abundance. Calculation of 
CRCs is most useful for species which uniquely trace cer- 
tain physics in a wide variety of astrophysical objects. 

As a second criterion, the possibility to use old results or 
results from other molecules is considered. Although direct 
calculations are always preferred over approximate scalings, 
the impact of new calculations is somewhat less if scaling 
may be expected to yield fairly accurate approximations. 
The most common examples of such scalings are isotopo- 
logues, where CRC scaling will only work if the symme- 
try of the molecule is unchanged (e.g., HCO+— >DCO^ will 
work but H2CO^HDCO will not); isomers, where scaling 
only works if the species have similar dipole moments (e.g., 
HCN^HNC works better than HCO+-»HOC+); and O^S 
substitutions (e.g., H2CO— >H2CS). Unlike in molecular physics, 
accuracy to order of magnitude is often enough in astro- 
physics, and relative rates are more important than absolute 
values. 

The astrophysically most important molecules without 
any known rate coefficients and without structurally similar 
species with known rates are: 



1 . CN is a key tracer of energetic radiation in protoplane- 
tary disks. Galactic PDRs and extragalactic systems es- 
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Fig. 1 Line ratios of HCO+ and HCN, calculated with RADEX in the optically thin limit 
density. Contours are spaced linearly and some contours are labeled for easy identification. 
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pecia lly when compared with HCN (IMeijerink & Spaans 
2OO5I). Rate coe fficients for electron collisions are known 



(iBlack & van D ishoeck 1991), but no values for H2 or 
He exist. Scaling from CO or CS fails because CN has a 
■^E ground state. Note that calculated rates may be scal- 
able to CO+, another tracer of high-energy radiation, re- 
cently detected in the active galaxy M 82 (iFuente et al 



centlv ' 
2006h . 



H2D+ and D2H+ are unique chemical and kinematic 
tracers of the centers of pre-stellar cores, whe re all CNO- 
bearing species are frozen out on dust grains jCaselli et al 



2003 



Van derTaket al.M2005l: iVastel et al.l 12004 ). The 



asymmetric light-weight structure prevents scaling from 
other molecules. 

NO is another readily observed important PDR and XDR 
tracer. Information from other ^FI molecules such as OH 
cannot be scaled because of the different Hund case. Ex- 
perimental cross sections have been measured by the 
group of Milliard Alexander in Maryland, and may be 
used as a starting point. The ra tes may be further sca led 
to NS, a possible shock tracer faatchell & Vitill2002l) . 
H3O+ is impo rtant as tracer of H2O a nd of the local 



ionization rate ( van der Tak et al. 20061). So far , scaled 



NH3 rates have been used (iPhillips et al.l 119921) which 
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are at best correct to order of magnitude. 



Finally it would be useful to develop approximations to 
derive rate coefficients for large molecules (27 atoms) for 
which detailed calculations would take prohibitive amounts 
of time. 
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